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CHAPTER  1 
INTRODUCTION 


IPROP  channel  tracking  simulations1  of  fast  rise  time  beams  in  conduc¬ 
tivity  channels  typically  predict  much  stronger  tracking  forces  than  those 
of  previously  published  computational  studies.2*3  Moreover,  the  forces 
are  predominantly  magnetic,  whereas  the  pivotal  analytical  model  In  use  had 
assumed,4  and  the  earlier  computational  studies  found,  only  electrostatic 
tracking  forces. 

Recent  Investigations5.6  suggest  that  the  larger  tracking  forces  in 
IPROP  are  in  part  due  to  differences  In  the  chemistry  model  used  and  the 
beam  profiles  assumed.  For  Instance,  IPROP  employs  a  temperature  dependent 
momentum  transfer  cross  section  which  Increases  by  as  much  as  a  factor  of 
three  as  the  channel  temperature  rises  due  to  the  presence  of  the  beam 
head.7  However,  significant  qualitative  and  quantitative  discrepancies 
between  the  tracking  results  of  IPROP  and  the  other  three-dimensional 
propagation  codes  persist  even  after  chemistry  and  beam  profiles  dif¬ 
ferences  are  taken  Into  account. 

These  remaining  discrepancies  seem  to  be  caused  by  the  different 
electromagnetic  field  models  used  in  the  codes.  IPROP8  solves  the  complete 
set  of  Maxwell's  equations  in  a  Galilean-transformed  frame  moving  with  the 
beam.  Alternatively,  IPROP  can  solve  the  Frozen  Field  Equations,  which  are 
obtained  from  Maxwell's  equations  by  assuming  the  fields  In  a  Galilean- 
transformed  frame  moving  with  the  beam  at  the  speed  of  light  are  Inde¬ 
pendent  of  time.  For  relativistic  beams  the  tracking  forces  are  essen¬ 
tially  the  same  In  both  cases.  RINGBEARER  II9  and  DYNASTY  II,10  on  the 
other  hand,  utilize  the  New  Field  Equations,11  an  approximation  to  the 
Frozen  Field  Equations  which  Includes  only  the  scalar  potential  and  axial 
component  of  the  vector  potential  with  no  displacement  current.  A  recent 
modification  to  DYNASTY  II  reintroduces  the  axial  displacement  current, 
which  often  increases  the  tracking  forces  seen  In  that  code.5 


To  clarify  the  impact  of  various  electromagnetic  field  approximations 
on  conductivity-channel  tracking,  we  have  developed  analytical  tracking 
models  for  the  Frozen  Field  Equations,  for  the  New  Field  Equations,  and  for 
the  New  Field  Equations  plus  the  axial  displacement  current.  As  in  the 
original  Electrostatic  Tracking  Model,4  the  channel  conductivity  is  assumed 
constant  in  space  and  time  within  the  channel,  which  has  a  sharp  edge,  and 
the  beam  is  contained  entirely  within  the  channel.  Although  somewhat  arti¬ 
ficial,  these  assumptions  seem  to  capture  much  of  the  essential  tracking 
phenomena  while  making  the  calculations  analytically  tractable.  Extending 
the  models  to  accommodate  smoothly  varying  radial  conductivity  profiles  is 
desi rable. 

The  Frozen  Field  Tracking  Model  reproduces  well  the  results  of  IPROP 
when  the  code  is  run  under  the  same  conditions.1  In  particular,  the  model 
correctly  predicts  strong  magnetic  tracking  and  electric  detracking,  with 
the  former  dominant,  very  early  in  the  beam  pulse.  As  a  certain  quantity 
X^,  to  be  defined  later,  goes  to  zero,  the  model  equations  reduce  to  those 
of  the  Electrostatic  Tracking  Model.  Even  for  moderate  the  total 
tracking  force  from  the  Frozen  Field  Tracking  Model  typically  agrees  with 
the  force  from  the  Electrostatic  Tracking  Model  to  within  a  factor  of  two, 
although,  of  course,  there  is  no  agreement  at  all  between  the  electric  and 
magnetic  components  of  the  forces  in  the  two  models.  Approximate  agreement 
between  the  total  forces  arises  from  the  near  cancellation  of  the  electric 
and  magnetic  forces  In  Frozen  Field  Tracking  Model.  Whether  the  agreement 
persists  for  more  realistic  channel  conditions  is  uncertain. 

Quite  surprisingly,  the  tracking  model  based  on  the  New  Field 
Equations  with  an  axial  displacement  current  gives  electric  and  magnetic 
tracking  force  equations  formally  identical  to  those  of  the  Frozen  Field 
Tracking  Model  but  with  X  different  by  about  a  factor  of  two.  The  tracking 
model  for  the  New  Field  Equations  without  an  axial  displacement  current,  on 
the  other  hand,  turns  out  to  be  the  same  as  the  Electrostatic  Tracking 
Model . 
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This  report  describes  in  detail  the  derivation  of  the  Frozen  Field 
Tracking  Model  and  presents  sample  predictions.  In  addition,  it  outlines 
corresponding  derivations  for  the  New  Field  Equations  with  and  without  an 
axial  displacement  current.  The  derivations  are  straightforward  but 
lengthy.  Source-free  dipole  fields  are  determined  in  the  regions  inside 
and  outside  the  channel.  The  beam  fields  then  are  derived  without  con¬ 
sideration  of  the  channel  boundary.  Matching  the  three  fields  at  the 
channel  boundary  determines  coefficients  appearing  in  the  source-free  solu¬ 
tions,  and  the  interior  source-free  solution  integrated  over  the  beam  cross 
section  gives  the  dipole  tracking  force.  Because  the  derivation  involves 
only  a  two-term  expansion  in  azimuthal  angle,  it  is  strictly  valid  only  for 
small  displacements  of  the  beam  from  the  channel  axis.  However,  IPROP 
simulations  suggest  that  the  models  are  accurate  for  larger  displacements, 
provided,  of  course,  that  the  beam  remains  in  the  channel. 

The  material  is  organized  as  follows.  In  Chapter  2  the  Electrostatic 
Tracking  Model  Is  rederived  using  the  procedure  just  outlined  in  order  to 
illustrate  this  procedure  in  a  simple  case  and  to  provide  a  basis  for  sub¬ 
sequent  comparisons  with  the  other  models.  Chapter  3  contains  the  main 
results  of  the  report.  The  Frozen  Field  Equations  are  collected  and  cast 
as  a  coupled  pair  of  two-dimensional  diffusion  equations.  For  analytical 
tractability,  the  diffusion  terms  then  are  dropped,  and  the  equations 
solved  to  obtain  the  electric,  magnetic,  and  total  tracking  forces. 

Various  limits  are  considered,  including  the  Electrostatic  Tracking  Model. 
This  same  approach  Is  applied  to  the  New  Field  Equations  in  Chapter  4. 

As  already  noted,  it  leads  to  the  Frozen  Field  Tracking  Model  with 
redefined  when  the  axial  displacement  current  is  retained  and  to  the  Elec¬ 
trostatic  Tracking  Model  when  it  Is  not.  Neglecting  the  diffusion  terms  is 
justified  in  Chapter  5.  The  diffusion  equations  are  Laplace-transformed  in 
the  axial  direction  and  solved  by  a  procedure  similar,  but  not  identical, 
to  that  used  in  earlier  sections.  A  low  frequency  expansion  then  recovers 
the  transforms  of  the  earlier  results  and  provides  an  estimate  of  their 
range  of  validity.  The  derivations  in  Chapters  2  and  3  also  assume  that 
the  radius  of  the  experimental  chamber  is  large  compared  to  the  radius  of 


the  channel.  The  effects  of  a  smaller  chamber  are  estimated  in  Chapter  6. 
Conclusions,  including  a  partial  physical  explanation  of  the  predicted 
tracking  behavior,  are  provided  in  Chapter  7. 

The  electron  charge  and  mass,  and  the  speed  of  light  are  set  equal  to 
unity  throughout  the  report.  A  factor  of  4w  is  absorbed  into  the  conduc¬ 
tivity  <3  and  the  beam  density. 


CHAPTER  2 

ELECTROSTATIC  TRACKING  MODEL 


Lee's  Electrostatic  Tracking  Model4  considers  a  rigid,  axisymmetric 
electron  beam  propagating  in  a  sharp  edged,  axisymmetric,  constant  conduc¬ 
tivity  channel  of  radius  a.  The  beam  is  entirely  within  the  channel  and 
offset  slightly  from  its  axis  by  a  distance  e,  as  depicted  in  Figure  1. 

(The  drift  tube,  of  radius  j),  plays  no  explicit  role  in  this  chapter.) 

Computations  are  performed  in  the  beam  frame,  assumed  to  move  at  the 
speed  of  light.  The  electrostatic  potential  is  determined  from 

vj2*  *  -0  +  pc  (1) 

p  and  p  are  the  beam  and  channel  plasma  charge  densities,  the  latter  given 
0 

by" 

It  bc  +  7l  •  *  o  (2) 

o  is  the  channel  conductivity,  x  represents  axial  distance  measured  back 
from  the  head  of  the  beam  (i.e.,  t  *  t  -  z),  and  r_  and  e  have  their  usual 
meanings  in  cylindrical  coordinates.  V£  Is  the  transverse  (i.e.,  £  and  0) 
gradient  operator. 

Combining  Equations  1  and  2  yields 

It  712*  +  |t  j  (3) 

where  p  has  been  replaced  by  J,  the  beam  axial  current  density,  for  nota- 
tional  consistency  with  subsequent  chapters.  (The  two  quantities  are  equal 
for  the  rigid,  ultrarelativistic  beam.)  Both  inside  and  outside  the 
channel,  where  o  is  constant.  Equation  3  can  be  integrated  once  in  x. 
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FIGURE  1.  ELECTRON  BEAM  IN  SHARP  EDGED  CONDUCTIVITY  CHANNEL  OF  RADIUS  a  AND 
DRIFT  TUBE  RADIUS  OF  b.  THE  BEAM  CENTROID  IS  DISPLACED  FROM  THE 
CHANNEL  AND  DRIFT  TUBl  AXIS  BY  A  SMALL  DISTANCE  e. 


J  dT‘ 


(4) 


% 


Vj>  -  -  /  e-a(T-,') 
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The  radial  derivative  of  ♦  may  be  discontinuous  at  the  channel  boundary  due 
to  surface  charges.  From  Equation  3, 


i_i_  (e) 
3t  3 r  v 


(5) 


The  potential  itself  is  continuous. 


♦ 


(e) 


(6) 


Equations  4-6  are  now  solved  for  the  force  acting  on  the  beam.  The 
dipole  (e.g.,  varying  as  cos8)  component  of  *  inside  the  channel. 


a  r  +  S 


(7) 


consists  of  two  terms.  is  the  dipole  potential  of  the  beam  in  a  uniform 
conductivity  background  without  a  boundary.  Thus,  it  contributes  nothing 
to  the  tracking  force.  Outside  the  beam. 


5.2l  JV‘ <",>1irI<lt. 
0 


(8) 


I ( t )  is  the  beam  current,  the  integral  of  J  over  the  beam  cross  section. 
The  second  term,  ar/a,  is  generated  by  surface  charges  at  the  channel 
boundary.  Outside  the  channel  the  dipole  potential  is  simply. 


Like  a,  y  is  an  unknown  function  of  t  to  be  determined. 


(9) 
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Inserting  Equations  7  and  9  Into  the  boundary  conditions.  Equations  5 
and  6,  provides  a  differential  equation  for  a. 

an) 

Its  solution  is 

a  *  — §  i  .  e-a(,-,’>)  ,dt,  (n) 

a  0 

The  corresponding  dipole  electric  fields  E_  and  Ea  are  give  by  -otcose  and 
-asine,  respectively.  Averaged  over  the  beam  cross  section  and  added,  they 
give  the  tracking  force,  -ol.  This  is  indeed  the  result  derived  by  a  dif¬ 
ferent  route  in  Reference  4.  (Actually,  Lee's  expression  is  slightly  more 
general  in  that  it  allows  o  to  vary  with  t.)  The  negative  sign  indicates 
that  the  force  is  inward.  Physically,  plasma  electrons  repelled  by  the 
beam  first  accumulate  preferentially  at  the  channel  edge  closest  to  the 
beam  at  the  rate  a  and  then  redistribute  themselves  uniformly  around  the 
channel  edge  at  the  rate  a/2 .  Tracking  persists  in  this  model  as  lony  as 
3I/3t  >  0. 


8 


•J  S 

*  '•  '  ■■  '■  1  •  •  *.;•••! 


CHAPTER  3 

FROZEN  FIELD  TRACKING  MODEL 


The  Frozen  Field  Equations  are  a  simplified  form  of  Maxwell's  equa¬ 
tions  obtained  by  performing  a  Galilean  transformation  to  a  frame  moving  at 
the  speed  of  light  and  then  dropping  time  derivatives. 


feEz-fl?rB9-FleBr-»E2-Jz 


1_R,  I  !_  rF  4.  II-  F 

3t  Bz  *  r  3r  E0  r  90  Er 


3t  (Er  B0^  ar  Ej 


—  (E  -  b  )  3  —  —  B  -oE  -  J 
3t  1  r  0J  r  30  z  r  r 


3t  (E0  +  V  ’  *  r  99  Ez 


It  (Ee  +  Br)  *  -  Ir  Bz  -  *E0  - 


(12) 

(13) 

(14) 

(15) 
(W 
(17) 


These  equations  must  be  used  with  care,  since  their  derivation  Involves  the 
Implicit  assumption  that  sources  move  at  the  speed  of  light.  Nonphysical 
behavior  sometimes  arises  near  Infinite  conductivity  surfaces,  and  electro¬ 
magnetic  shock  waves  are  possible.  Nonetheless,  the  Frozen  Field  Equations 
have  been  used  successfully  In  many  particle  beam  propagation  calculations. 
As  noted  in  the  Introduction,  IPROP  channel  tracking  simulations  show 
negligible  change  (except  for  beams  passing  through  metal  foils)  when 
Maxwell's  equations  are  replaced  by  the  Frozen  Field  Equations. 


For  a  *  0,  It  Is  convenient  to  combine  Equations  12-17  Into  a  pair  of 
second  order  equations  for  Ez  and  Bz,  analogous  to  the  usual  equations 
describing  electromagnetic  fields  in  waveguides. 


S 


*•‘1 


- -F To i If] BI  -^J2  <18> 

1_  n  a  \l  i_  -  Ii_  .  1  jLilil  O 

3t  z  I  r  Jr  o  ir  r  36  a  r  30  J  z 

-j"ii_r— -- — I  i_  I  A_1  £  (19) 

[r  Jr  j  r  39  r  30  a  3r J  z 

Transverse  beam  currents,  J„  and  J  ,  have  been  dropped,  consistent  with  the 

r  o 

assumptions  presented  in  Chapters  1  and  2.  Once  the  longitudinal  fields 
are  determined,  the  transverse  field  are  obtained  from 

Er-  i  (If  Ez  *  F  le  8l)  <20> 

E»*  -lAb*z-7h*z)  <21> 

It  Br  ■  If  -  (ltT?i)rliE«  <22> 

fr  -  (l  +  TF  i)  If  E*  *  I?  i  f  I*  8z  <23> 

Equations  18  and  19  simplify  greatly  for  constant  conductivity.  If  a 
Is  independent  of  r_  and  6,  the  equations  decouple.  If,  In  addition,  a  is 
Independent  of  t.  Equation  18  can  be  Integrated  once. 

?  1  x 

_  a  l  a  a  l  a4  r  ,  „  f  -o( t-t1  )  a  ,  ...  (0.s 

FTF  r1r  '  7  ^zj  Ei  *  J  e  W  J  dT  (24> 


i'-,  r  *  r  n  r 


Both  are  two-dimensional  diffusion  equations.  Note  the  similarity  between 
Equation  24  and  Equation  4. 

Equations  24  and  25  can  be  solved  exactly  Inside  the  channel  using 

Green's  functions  or  Laplace  transforms.  However,  the  formal  solution 

which  results  Involves  integrals  with  complicated  kernels  and  is  not 

particularly  enlightening.  Instead,  we  obtain  an  approximate  solution  by 

dropping  the  diffusion  terms  o3/3t  and  treating  the  simplified  equations  by 

the  method  used  in  the  preceding  Chapter.  Omitting  the  diffusion  terms  in 

Equations  24  and  25  is  reasonable,  if  fields  change  in  t  on  a  scale  slow 
2 

compared  to  aa  .  The  range  of  validity  of  this  approximation  is  estimated 
in  Chapter  5,  where  the  Laplace  transformation  solution  is  obtained  and 
then  expanded  in  powers  of  the  transform  variable.  Comparing  the  approxi¬ 
mate  analytical  solution  with  IPROP  computations  presented  later  also 
supports  dropping  the  diffusion  terms. 

With  the  ct3/3t  omitted.  Equations  24  and  25  become 

V1  Ez  “  <26> 

0 

4bzs°  (27) 

Their  dipole  solutions  are 

4°  *  It  or  “  aS  (28) 

B^-f^Sr  (29) 

E,  is  taken  to  vary  as  cose.  In  which  case  E_  and  B.  also  vary  as  cose,  and 
z  r  o 

B  ,  Ea  and  B„  vary  as  sine.  Derivatives  of  the  x-dependent  unknown  func- 

Z  u  P 

tions  a  and  8,  rather  than  the  functions  themselves,  are  Introduced  to 
simplify  notation  later  in  the  calculation.  S  was  introduced  in  Chapter  2 
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and  is  expressed  by  Equation  8  for  that  portion  of  the  channel  outside  the 
beam. 


Equations  18-23  cannot  be  used  outside  the  channel,  because  their 
derivation  involves  division  by  a.  Simply  multiplying  those  equations  by 
a  before  setting  it  to  zero  leaves  an  incomplete  system  of  equations.  In 
essence,  the  Frozen  Field  Equations  are  qualitatively  different  in  vacuum 
and  in  resistive  plasma. 


Rearranging  Equations  12-17  to  evaluate  the  fields  outside  the  channel 
proceeds  as  follows.  As  before,  the  longitudinal  field  components  are 
easily  isolated. 


v|  Ez  =  vT  8. 


d 


(30) 


However,  Ez  and  Bz  no  longer  are  independent.  From  Equations  14  and  15, 


+  1 b  -  0 

3r  z  r  30  z 


(31) 


In  consequence,  the  transverse  field  components  are  only  partly  determined 


by  Ez  and  Bz. 


i.(j  *r  *  J  !«  E.)  ’  '  7  E2 


ii-rl  E 
r  ir  )r  ! 


(32) 


(33) 


(34) 


(35) 
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.  VYl  'Ifu 


mgit 


From  Equations  30  and  31,  the  dipole  longitudinal  fields  are 


(36) 


2  2 

if  a  /b  c  1.  Of  the  dipole  transverse  fields,  only  the  combination 
Ea  -  B0  is  needed  to  derive  the  tracking  force  In  the  large  radius  drift 

9  9 

tube  limit.  It  is  described  by 


3t  dr  (Ee  *  Be) 


(37) 


which  is  the  sum  of  Equations  32  and  33  evaluated  for  the  azimuthal  field 
dependences  given  following  Equation  29. 


E(e)  _  Bie)  =*  2a2  - 


-  v 


(38) 


Note  that  the  drift  tube  radius  J>  appears  logarithmically  in  the  first  term 
of  Equation  38.  That  term  is  not  necessarily  small  and  cannot  be  dropped. 
The  vacuum  fields  are  derived  for  arbitrary  b/a  in  Chapter  6. 


Three  matching  conditions  at  the  channel  boundary  are  necessary  to 
determine  the  independent  quantities  a,  8  and  y. 


(39) 


B<e>  «  ^ 


(AO) 


I  !_  1  !_  r  (e(1) 

a  3t  r  3r  V  z 


i_  .(I) 
3r  Lz 


(41) 


Equation  41  follows  from  continuity  of  Ea  -  8tt,  and  from  Equations  21  and 

0  0 

23. 
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The  jump  conditions  Equations  39-41  can,  of  course,  also  be  treated  as 
boundary  conditions  on  the  interior  solutions  by  eliminating  the  exterior 
solutions. 


r(i)  ,  B(i) 


2a2  tn  b/a  i-*  E^ 1 J 
3t 


»  i-il_  r  +  8{i)^  +  a  li-  rE 

a  3:  p  3r  r  \Lz  Hz  /  +  a  r  Tr  rt 


(i) 

z 


(42) 

(43) 


The  tracking  force  is  now  obtained.  Substituting  Ez  and  Bz  from 
Equations  28,  29,  and  36  into  Equations  39  and  40  gives  B  *  y  and 


(44) 


Equations  28,  29,  38,  41  and  44  then  yield  a  differential  equation  for  a. 


where 

X  s  o2a2  tn  b/a  (46) 

The  solution  to  Equation  45  is 


'7  j  {*'<’<T",')/X  [cosh^3.(,  -  ,•) 


+ - - —  slnh  —  —  o( t  -  t') 

x 


g-o(  t-t'  ) 


(47) 


The  net  angle-averaged  dipole  field  acting  on  the  beam  Is  1/2(E  -  Ba  - 

r  0 

Efl  -  Br).  This  quantity  Is  expressed  in  terms  of  Ez  by  using  Equations  14 
and  16. 


B9  -  Ee  -  M 


-  —  r  E 
r  3r  z 


(48) 
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As  In  Chapter  2,  the  term  proportional  to  _S  in  Ez  does  not  contribute  to 
tracking.  Hence,  the  Frozen  Field  Tracking  Model  force  is  -ol. 


For  subsequent  comparison  with  IPROP  data,  the  separate  electric  and 
magnetic  tracking  force  components  are  useful.  The  transverse  fields 
inside  the  channel  are  obtained  from  Equations  20-23. 


r  0  a  3t 


Br  =  89  =  ‘  7  77 6 


(49) 


(50) 


Inserting  a  and  8  from  Equations  44  and  45  then  gives 


|C'r  /  I''’1"' 


-t')/x  2  /m 


sinh 


a(t  -  t'  ) 


-«(t-t')1  3 

|  TP 


I  dt' 


(51) 


P) 


- l—  sinh  ^37  a(T  .x')  Xl  dt'  (52) 

✓n  x  J 3T 

Multiplied  by  -I,  these  are  the  electric  and  magnetic  forces,  respectively. 
It  is  evident  from  Equations  51  and  52  that  the  electric  force  is  detrack¬ 
ing,  and  the  magnetic  force  tracking,  at  early  times. 

It  is  clear  from  Equations  10  and  45  that  the  Frozen  Field  Model 
tracking  force  must  reduce  to  the  Electrostatic  Model  tracking  force  as 
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X  ♦  0.  The  magnetic  component.  Equation  52,  vanishes,  while  the  electric 
component.  Equation  51,  approaches  the  electrostatic  force.  Equation  11. 
This  convergence  Is  not  uniform,  however.  Tracking  at  the  very  front  of 
the  bean  Is  modified  by  magnetic  effects  for  a  distance  in  t  of  order  X/o 
no  matter  how  small  X  becomes.  Of  course,  strong  modifications  over  very 
short  distances  have  little  practical  Importance. 

Although  perhaps  not  Immediately  obvious,  the  tracking  force  expres¬ 
sions  are  continuous  through  X  *  1.  At  X  *  1,  Equation  47  becomes 

.  j  e'<’(T-T')  a(T  -  t')  I  dt'  (53) 

1  0 

The  analytic  extension  of  Equation  47  above  X  *  1  Is  obtained  by  replacing 
(1  -  X)1^2  by  (X  -  l)^2  and  the  hyperbolic  functions  by  the  corresponding 
trignometrlc  functions.  Thus,  the  tracking  force  can  oscillate  in  ampli¬ 
tude  for  X  >  1.  Eventually,  dropping  the  diffusion  terms  from  Equations  24 
and  25  becomes  Inaccurate  as  X^  is  increased. 

The  variation  of  the  predicted  tracking  force  with  X^  for  a  10  kA, 

15  cm  rise  length  beam  offset  0.25  cm  In  a  1  cm  radius,  1  cm-1  constant 
conductivity  channel  is  Illustrated  In  Figure  2.  The  force  is  nearly  inde¬ 
pendent  of  ,X  and  equal  to  the  Electrostatic  Tracking  Model  force  for  X  <  2. 
A  stronger  but  oscillatory  force  arises  at  larger  X.  Note  that  the  25-50  g 
peak  tracking  forces  should  not  be  expected  for  more  realistic  beam  and 
channel  profiles. 

The  similarity  between  tracking  force  predictions  from  the  Electro¬ 
static  and  Frozen  Field  models  suggested  In  Figure  2  does  not  carry  over  to 
their  electric  and  magnetic  components.  Figure  3  depicts  the  tracking 
force  components  from  Equations  51  and  52  for  the  beam  and  channel  of 
Figure  2  but  with  o  ■  0.5  cm"1;  X  *  1.  For  a  distance  into  the  pulse  of 
order  a  few  times  X/o,  the  magnetic  component  clearly  dominates.  In  fact, 
the  electric  component  is  detracking  for  a  short  distance.  The  early-time 
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TRACKING  FORCE  (GAUSS) 


FIGURE  2. 


R-899 


THE  FROZEN  FIELD  TRACKING  MODEL  NORMALIZED  FORCE  (-a  FROM  EQ.  47) 
ON  A  10  kA,  15  an  RISE  LENGTH,  0.5  cm  RADIUS  BEAM  OFFSET  e  = 

0.25  cm  IN  AN  a  -  1  cm,  a  *  I  cm"!  CONSTANT  CONDUCTIVITY  CHANNEL 
FOR  X  ■  0,  1,  2,  AND  5.  THE  X  •  0  LIMIT  CORRESPONDS  TO  THE 
ELECTROSTATIC  TRACKING  MODEL  RESULT. 


TRACKINC  FORCE  (GAUSS) 


FIGURE  3.  THE  ELECTRIC  AND  MAGNETIC  COMPONENTS  OF  THE  NORMALIZED  TRACKING 
FORCE  ON  A  10  kA,  15  cm  RISE  LENGTH,  0.5  cm  RADIUS  BEAM  OFFSET 
e  -  0.25  cm  IN  AN  a  =  1  cm,  a  *  0.5  cm-1  CONSTANT  CONDUCTIVITY 
CHANNEL  FOR  X  »  1.  THE  FORCES  ARE  OBTAINED  BOTH  FROM  EQUATIONS 
51  AND  52  AND  FROM  A  COMPARABLE  IPROP  SIMULATION. 


magnetic  tracking  appears  to  be  due  both  to  lateral  plasma  currents  and  to 
the  axial  displacement  current.  Inductive  effects  cause  the  electric 
tracking  reversal  at  the  beam  front.  The  remarkable  cancellation  between 
the  components  which  allows  the  total  tracking  force  to  so  nearly  equal  the 
purely  electrostatic  result  from  Chapter  2  may  not  persist  for  more 
realistic  channel  profiles. 

Also  plotted  in  Figure  3  are  the  tracking  forces  from  a  comparable 
IPROP  simulation.  Agreement  is  very  good  overall,  supporting  the  validity 
of  our  tracking  model.  Radial  profiles  of  the  various  fields  in  IPROP 
agree  well  with  the  predictions  of  Equations  20-23.  IPROP  preliminary 
tracking  results  are  presented  in  some  detail  in  Reference  1. 


CHAPTER  4 

CHANNEL  TRACKING  WITH  THE  NEW  FIELD  EQUATIONS 


Deriving  the  tracking  force  implied  by  the  New  Field  Equations11  is 
desirable,  because  this  approximation  to  the  Frozen  Field  Equations  is  used 
in  many  propagation  codes.  The  New  Field  Equations  can  be  written  in  many 
equivalent  forms,  including3 

~  v|(Az  -  ♦)  *  Vj_  •  aVj_*  (54) 

2 

tAz  -  ♦)  +  «  It  (A2  -  ♦)  -  4  A2  -  J  (55) 

3  T 

The  axial  displacement  current,  sometimes  used  in  DYNASTY  II,10  is  added  to 
Equation  55.  From  time  to  time  during  the  derivation  which  follows,  we 
note  the  effects  of  dropping  this  term.  In  this  way  the  tracking  forces 
for  the  New  Field  Equations  with  and  without  the  axial  displacement  current 
are  obtain  simultaneously. 

Inside  the  channel,  where  o  is  constant,  Equations  54  and  55  can  be 
manipulated  to  provide  an  equation  for  Az  -  ♦  alone. 

Mir  *  •)(*!  -  ■  (l  A  (A2  -  ♦)  ♦  Jz  m 

Either  Az  or  +  can  be  taken  as  the  second  dependent  variable.  We  choose 

v 

iAz-  (»  *?fr)  4  CAZ  -  ♦)  (57) 

As  in  the  preceding  derivations.  Equation  56  Is  simplified  by  Integrating 
once  in  t. 
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('k- 4)  (*.-♦)■«  J 


(58) 


Note  that  the  first  term  in  Equation  58  has  a  more  complicated  form  in  the 
absence  of  the  axial  displacement  current.  However,  as  before,  we  drop 
this  diffusion  term.  In  other  words,  the  axial  displacement  current  plays 
no  direct  role  inside  the  channel. 


T 

v|  (Az  -  ♦)  -  -O  J*  e'a(T‘T,)  J  dt* 


Equations  59  and  57  have  dipole  solutions 


(59) 


-♦(1) 


=  -ar  +o  I  S  dt' 


T 

J 

0 


411  ■  (‘  *14?)  (411  -  ♦(,))* 8r 


(60) 


(61) 


In  the  vacuum  region  outside  the  channel.  Equations  54  and  55  reduce 


to 


v|  (Az  -  ♦)  -  0 

7i  Az  “  7T  (Az  ‘ 

O  T 


Their  dipole  solutions  are 


=  5  y  in  b/r 

3t 


(62) 


(63) 


(64) 


(65) 
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As  in  Chapters  2  and  3,  small  terms  scaling  as  (a/b)fc  are  omitted.  The 
second  term  in  Equation  65  arises  from  the  axial  displacement  current. 


Four  matching  conditions  at  the  channel  boundary  are  needed  to 
evaluate  a,  8,  y»  and  5.  Az  and  $  are  continuous  there.  Integrating 
Equations  54  and  55  across  the  boundary  gives,  in  addition. 


—  A(e) 
3r  *z 


i_  A(1) 
3r  Az 


!_  i_  /.(e)  _  .(e)\  s  3_  3_  /.(i)  _  (i)\  3 

3t  3r  \Az  ♦  )  3t  3r  \Az  ♦  )  alr 


,0) 


(66) 

(67) 


Combining  these  with  Equations  64  and  65  to  eliminate  the  exterior  fields 
yields 


If  desired,  A^  also  can  be  eliminated,  leaving  a  single  boundary 
condition  at  r  *  a. 


(68) 

(69) 


•(“  b/a  - 1/2)  77  (41’  -  ♦(1))  +  i  (Sr ♦  'Hfe r  (Azn  -  *(1>) 

The  second-derivative  terms  In  Equations  68  and  70  are  omitted  In  the  case 
of  no  axial  displacement  current. 


Inserting  Equations  60  and  61  Into  Equations  68  and  69,  or  Equation  60 
Into  Equation  70  yields  a  differential  equation  for  o  Identical  to  Equation 
45  for  the  Frozen  Field  Tracking  Model  but  with  X  redefined  as 
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X  ■  |  c2a2(w  b/a  -  1/2) 


(71) 


With  no  axial  displacement  current,  X  ■  0.  The  tracking  fields  are 
obtained  by  taking  the  gradient  of  -ar«cos§.  When  averaged  In  angle  and 
summed,  they  give  a  tracking  force  -ml,  as  before.  For  completeness,  we 
note  that  the  electric  and  magnetic  components  of  the  tracking  force  are 
given  by  Equations  51  and  52  with  £  fro &  Equation  71. 

In  summary,  the  New  Field  Equations  give  rise  to  a  tracking  force 
equal  to  that  of  the  Frozen  Field  Tracking  Model  (with  X  reduced  by 
approximately  one-half)  or  that  of  the  Electrostatic  Tracking  Model, 
depending  on  whether  the  axial  displacement  current  Is  Included.  This 
progression  X  ♦  1/2X  ♦  OX  Is  associated  with  dropping  first  3Bz/3t  and  then 
3Ez/3t  from  the  vacuum  field  equations  and  probably  Is  due  to  the  reduced 
field  energy  there.  The  plasma  currents  both  drive  and  are.  In  part, 
driven  by  this  field  energy,  which  Is  proportional  to  tn(b/a). 
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CHAPTER  5 

DIFFUSION  CORRECTIONS 

i 

! 

Diffusion  terms  were  ignored  in  previous  sections  in  order  to  obtain 
analytically  tractable  results.  Here,  we  evaluate  the  Laplace  transforms 
of  the  tracking  forces  with  the  diffusion  terms  retained  in  order  to  esti¬ 
mate  their  importance.  The  transformed  tracking  forces  can  be  derived  by 
paralleling  the  derivations  in  Chapters  3  and  4,  starting  from  the  trans¬ 
forms  of  Equations  24  and  25  for  the  Frozen  Field  model  and  Equation  58  for 
the  New  Field  models.  Indeed,  we  have  done  so  and  obtained  the  results 
given  below.  Here,  however,  we  outline  an  alternative,  more  formal 
approach. 


Throughout  this  section,  all  quantities  are  taken  to  be  Laplace- 
transformed  with  respect  to  t.  The  transform  variable  is  s. 


In  coordinates  centered  on  the  channel  axis,  the  net  tracking  force 
consists  of  a  monopole  force  acting  on  the  beam  dipole  current  and  a  dipole 
force  acting  on  the  beam  monopole  current.  The  dipole  current  Is  given  by 
-e  3J/3r  for  a  rigidly  displaced  beam.  Thus,  the  total  force,  averaged 
over  angle  and  normalized  to  the  beam  current.  Is 

rl|  [F0  (-'#)♦  F1JW  (72) 

The  forces  can,  in  turn,  be  expressed  in  terms  of  potentials. 


3r 


7IFr  51 


(73) 


(74) 


Inserting  Equations  73  and  74  Into  Equation  72  and  integrating  by  parts 
gives 


-I 


-1 


(75) 


^  7 


24 


The  integrand  of  Equation  75  would,  of  course,  vanish  identically  were  it 
not  for  the  channel  edge. 


Now  let  us  suppose  that  the  potentials  satisfy 

,(m)  o 


/i  iL  c  i_  .  k2  -4V 

\r  3r  r  3 r  J Si 


*  -J' 


a  +  s 


(76) 


with  boundary  conditions 
c0  +  n0  3r  c0  "  0 

h  +  nl7l?f  51*° 

at  the  channel  wall,  r  -  a.  These  equations  are  solved  conveniently  in 
terms  of  Green's  functions. 


(77) 


(78) 


4 m 


-Trhr  r'dr' 


(79) 


r  >  r' 
r  <  r‘ 


(80) 


Yq(  <a)  -  nQKY^K*) 
*0  *  Yo(,er)  ’  J0(Kr)  JQ(iea)  -  r^TOT 

Yj(iea)  +  t,i,cYq(  <a) 

♦l  "  -  Jl(,er)  J^wa)  ♦  n^J pCTiT 


(81) 


(82) 


J  and  Y  are  Bessel  functions  of  the  first  and  second  kinds.12 
m  m 

The  Integrations  In  Equations  75  and  79  are  straightforward,  if  the 
•  beam  profile  is  square,  so  that  3J/3r  *  -0  6( r  -  a^). 


25 


(83) 


it  a 
7  a  +  s 


) 

1  1/2  Va”  *0(,cr) 


r^ab 


J, («aK) 

it  a  T  1  d  .  /  _\ 

7  7TT  1  Ke  T/TTT  ¥‘r) 


ria  b 


(84) 


The  tracking  force  is 


2eg  I  /JjE^bA  _ t1  *  "O'1! 

’  *  5  a2 \1/2  ‘V  [Jjlrt)  +  ^KJjtrt) 


k)  1/2  Ki 


)3CJ0(-ca)  -  HgicJ j(  <a)] 


(85) 


The  Bessel  function  Wronskian  formula12  was  used  in  obtaining  Equation  85. 


It  can  be  shown  that  ng  *  0  for  the  three  tracking  models  Investigated 

in  this  report.  Physically,  this  corresponds  to  the  vanishing  of  the 

axisymmetric  portion  of  Ez  outside  the  channel  for  ultrarelativistic  beams. 

Equation  85  can  be  simplified  notationally  by  absorbing  the  shape  factor 
2 

[J1(Kab)/(l/2Kab)]  into  I.  This  quantity  Is  equal  to  unity  at  long  trans¬ 
verse  wavelengths,  and  decreases  rapidly  for  l/2iea^  >  2.4.  With  these 
changes.  Equation  85  becomes 

2ea  I _ 1/2  xa 

"  o  +  s  fl2  JQ( Jx  +  njxOg)  K  1 

From  here  on,  the  Bessel  function  arguments  are  understood  to  be  ica. 

For  the  Frozen  Field  Equations  of  Chapter  3,  -E  /s  can  be  identified 

2  z 

as  the  potential  and  -as  as  <  .  See  Equations  24  and  48.  Since  Bz  is 
proportional  to  Jm(«a)  Inside  the  channel.  Equations  42  and  43  can  be 
combined  to  yield  the  boundary  condition 

(2Xs2  +  os  Jq/Jj)  Ez  +  oa(o  Ml^rE^O  (87) 

from  which  can  be  Identified.  Consequently,  the  Laplace-transformed 
tracking  force  is 
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Similarly,  A  -  4  is  the  potential  for  the  New  Field  Equations  of 
Z  2 

Chapter  4.  The  transverse  wavenumber  <  is  -os  if  the  axial  displacement 
current  is  retained,  and  -os/(o  +  s)  otherwise.  See  Equation  58  and  the 
associated  text.  The  Lapl ace-transform  of  Equation  70  provides  the 
boundary  condition.  The  resulting  tracking  force  is 

2e  s  l  _ (Xs  ♦  o)q _ 

a2  a  +  s  JQ[(Xs2  +  os)J1  /  1/2 tea  +  a(a  +  s)JQ] 


Recall  that  X  Is  set  to  zero,  if  the  axial  displacement  current  Is  dropped. 


Setting  JQ  and  /  l/2iea  to  unity  In  either  Equation  88  or  89  yields 


2e  s 
,2  a  +  s 

fl 


(Xs  +  a) a 
Xs2  +  2os  +  a1 


(90) 


This  Is  the  Laplace  transform  of  Equation  47,  the  channel  tracking  force 
derived  with  the  diffusion  term  Ignored  Inside  the  channel,  as  can  be  seen 
most  readily  by  examining  Equations  8  and  45.  Thus,  dropping  the  diffusion 
term  Is  valid  if  js^J  «  4/aa^,  or  4o/(o2a^  +  1)  for  the  New  Field  Equations 
without  the  axial  displacement  current,  s^  is  a  zero  of  Xs2  +  2as  +  q2. 

Numerical  analysis  of  the  poles  and  residues  of  Equations  88  and  89  Indi- 

2 

cates  that  the  factor  4  can  be  Increased  to  2.4  . 


In  Chapters  3  and  4  and  again  here.  It  Is  seen  that  the  tracking 
forces  for  the  Frozen  Field  Equations  and  for  the  New  Field  Equations  with 
the  axial  displacement  current  are  formally  identical  when  the  diffusion 
term  Is  omitted,  differing  only  In  the  definitions  of  X.  The  complete 
force  expressions,  Equations  88  and  89,  are  not  formally  Identical,  how¬ 
ever.  Nonetheless,  their  lowest  poles  and  corresponding  residues  are 
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2  2 

nearly  equal  numerically,  at  least  for  o  a  <  10.  The  tracking  force  for 
the  New  Field  Equations  without  the  axial  displacement  current  is,  of 
course,  qualitatively  different  from  that  of  the  other  models,  unless 
X  «  1:  the  lowest  poles  are  a  complex  conjugate  pair  in  the  first  two 
cases  and  a  single  real  value  In  the  third. 

Further  numerical  Investigation  of  the  poles  and  residues  yields  a 
result  which  is  surprising  at  first  sight.  The  lowest  poles,  in  terms  of 
Re(s.j),  contribute  to  tracking,  or  in  the  case  of  complex  conjugate  pairs 
contribute  alternately  (in  t)  to  tracking  and  detracking,  but  tracking 
first.  The  lowest  pole  associated  with  detracking  comes  from  the  smallest 
zero  of  the  factor  JQ  in  the  denominators  of  Equations  88  and  89  and  so  has 
a  large  Re(s^).  (It  represents  beam  repulsion  from  the  channel  due  to  the 
monopole  return  current.)  This  seems  to  predict  that  beams  always  track 
conductivity  channels,  independent  of  the  magnitude  of  oa.  This  apparent 
discrepancy  with  simulation  results  probably  Is  resolved  as  follows.  The 
complex  conjugate  pair  which  occurs  for  X  *  0  has  a  very  large  Imaginary 
part  for  aa  very  large.  Thus,  the  tracking  force  oscillates  rapidly  in 
sign,  so  that  its  average  value  is  small.  In  addition,  the  corresponding 
residues  become  very  small.  When  X  *  0,  on  the  other  hand,  the  lowest 
tracking  pole  moves  very  close  to  the  first  detracking  pole  as  oa  becomes 
large,  and  their  effects  approximately  cancel.  The  behavior  of  higher 
poles  then  must  be  considered.  We  have  not  done  so,  because  It  does  not 
seem  worthwhile. 

Incidentally,  the  next  higher  tracking  pole  In  Equations  88  and  89 
2  2 

takes  on  the  value  3.8  /oa  for  oa  small  and  gradually  decreases  to 
2  2 

2.4  /oa  as  oa  becomes  very  large.  It  is  probably  this  pole  that  enters 
into  the  resistive  hose  instability. 
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CHAPTER  6 

CONDUCTING- WALL  CORRECTIONS 


The  conducting  drift  tube  surrounding  the  beam  channel  enters  the 

preceding  calculations  only  through  the  logarithm  appearing  in  X.  Other 

2 

terms,  of  order  (a/b)  or  smaller,  are  ignored.  Now,  we  repeat  the 
analysis  of  Chapter  3  with  these  terms  included  to  demonstrate  that  they 
are,  indeed,  unimportant  for  large  drift  tubes.  The  impact  of  a  small 
drift  tube  on  channel  tracking  also  is  determined. 


General  solutions  to  Equations  30  and  31  for  Ez  and  Bz  in  vacuum 
between  £  and  £  with  Ez(b)  *  0  are 


E<«)  ,  i.  V  a(l  .  4^ 

*  *  \r 

B(e) 

r  \r  *  7  *  / 


(91) 


(92) 


The  transverse  fields  are  derived  most  easily  by  taking  the  sums  and 
differences  of  Equations  32  and  33  and  of  Equations  34  and  35,  and  then 
integrating  once  each  In  t  and  r.  E,(b)  »  0  determines  one  Integration 
constant. 
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Matching  these  exterior  solutions  to  the  Interior  solutions,  already 
given  in  Equations  28  and  29,  is  accomplished  through  continuity  of  Ez 
(Eq.  39),  Bz  (Eq.  40),  E0  -  BQ  (Eq.  41),  and  E0  +  B0 


(e<«>  *  .*•»)  -  r  Ip  A  (Eo>  -  .o») 


♦  f  e«> 

3r  z 


(97) 


The  counterpart  of  Equation  97  was  not  required  in  Chapter  3  due  to  a 
fortuitous  cancellation.  Combining  all  these  equations  yields  after  some 
algebra 


X  3 
"  3l2 


2!^ 


a  3t 


2e 


(98) 


is  again  redefined. 

X  =  a2a2[in  b/a  +  1/4(1  -  a4/b4)] 


(99) 


The  primary  difference  between  this  formula  and  Equation  45,  which  was 
obtained  in  the  limit  of  large  b/a,  is  seen  to  be  the  additional  source 
term,  proportional  to  the  beam  current,  on  the  right  side  of  Equation  98. 
Other  terms  involving  J>  can.  In  fact,  be  eliminated  formally  by  renor¬ 
malizing  o,  X,  a,  and  t.  The  corrections  Involved  are,  in  any  event,  small 
2 

so  long  as  (b/a)  »  1.  Solving  Equation  98  yields  the  additional  contribu¬ 
tion  to  the  tracking  force. 


2c 

7 


f  g-a( t-t' )/X _ 1 

J  /n 


sinh  —  y  —  o( t  -  t')  crldt' 


(100) 


where  terms  of  order  (a/b)*  have  been  omitted.  For  X  «  1  this  becomes 
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(101) 


2|  j 
b  0 

2 

Both  expressions  reduce  to  2el/b  ,  the  usual  magnetic  image  force,  for 
oTr  »  2.  rr  is  the  current  rise  length.  The  corresponding  electric  image 
is  shielded  by  the  channel  conductivity  at  late  times.  Recall  that  these 
results  assume  oa  small  in  the  sense  of  Chapter  5.  The  magnetic  image 
force  is  reduced  by  return  currents  for  oa  large. 

2 

The  criterion  oTp  >  2  for  the  magnetic  wall  force  to  reduce  to  2el/b 
is  almost  always  satisfied  in  practice.  On  this  basis,  it  can  be  shown 

that  the  wall  force  on  the  beam  head  is  small  compared  to  the  tracking 

2 

force  when  (b/a)  »  10.  That  the  wall  force  on  the  beam  body  be  small  com- 

2 

pared  to  the  tracking  force  requires  that  (b/a)  be  much  greater  than 

several  times  or  ,  with  the  conductivity  evaluated  in  the  beam  head.  IPROP 
r  2 

simulations  suggest  (b/a)  »  500  for  realistic  beam  and  channel  radial  pro¬ 
files  with  a  -  1  cm,  a  -  1  cm-1,  and  -  100  cm.  It  is  desirable  that  the 
second,  more  restrictive,  condition  be  met  so  that  wall-force  centering  of 
the  beam  body  not  be  confused  with  tracking  of  the  beam  head. 
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CHAPTER  7 
CONCLUSIONS 


The  foregoing  calculations  provide  an  improved  model  of  conductivity 
channel  tracking.  The  derivation  is  based  on  the  complete  set  of  Frozen 
Field  Equations.  It  is  valid  for  oa  <  2.4  and  b/a  >  10,  although  b/a  »  20 
probably  is  necessary  for  unambiguous  experimental  confirmation.  The  model 
also  requires  the  beam  to  lie  entirely  within  a  sharp-edged  constant  con¬ 
ductivity  channel.  Note  surprisingly,  simulations  show  somewhat  weaker 
tracking  forces  for  more  realistic  conditions.  Qualitative  features  of  the 
model  nonetheless  persist. 

The  analytical  tracking  force  expression  contains  a  quantity  which 
measures  the  electromagnetic  tracking-field  energy  outside  the  channel. 

For  X  <  1  the  tracking  force  is  approximately  that  of  the  Electrostatic 
Tracking  Model.  For  X  >  1,  on  the  other  hand,  the  tracking  force  oscil¬ 
lates  in  time  like  an  underdamped  pendulum.  This  field  energy  can  be 
thought  of  as  inertia  of  the  tracking  force.  Dropping  Bz,  which  is  done  in 
the  New  Field  Equations  with  the  displacement  current  retained,  cuts  this 
inertia  in  half.  Dropping  the  displacement  current  as  well  eliminates  the 
inertia  entirely. 

^  diverges  logarithmically  for  large  drift  tube  radii.  This  seemingly 
nonphysical  result  is  a  consequence  of  the  ultrarelativistlc  approximation 
Implicit  in  the  Frozen  Field  Equations.  A  Lorentz  transformation  to  the 
beam  frame  shows  in  a  qualitative  sense  that  lb  should  be  replaced  by  yx^  in 
open  air.  rQ  is  the  axial  scale  length  over  which  the  tracking  force 
varies;  It  Is  smaller  than  xp.  Fields  tend  to  fall  off  spherically  rather 
than  cylindrlcally  at  radii  greater  than  ytq.  Values  of  tn(yx0/a)  larger 
than  about  15  are  unlikely  in  practice. 

A  distinctive  feature  of  the  Frozen  Field  Tracking  Model  is  the  role 
of  transverse  magnetic  fields.  Radial  and  azimuthal  dipole  plasma  currents 
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near  the  beam  head  drive  Bp  and  B0  via  Faraday's  law  provided  Bz  is  not  set 
to  zero.  In  addition,  the  dipole  axial  displacement  current  generates 
transverse  magnetic  fields  Irrespective  of  whether  Bz  =  0.  It  is  not  dif¬ 
ficult  to  see  that  Vx B  always  leads  to  a  centering  force  on  the  beam  at 
early  times,  when  the  lateral  plasma  currents  are  still  rising.  Moreover, 
the  force  can  be  quite  large.  That  the  transverse  electric  fields  should 
be  modified  inductively  so  that  they  largely  cancel  the  magnetic  force  for 
X  <  1,  leaving  a  net  tracking  force  not  much  different  from  that  of  the 
Electrostatic  Tracking  Model,  is  surprising.  Whether  this  cancellation 
persists  for  more  realistic  channel  profiles  is  unknown. 

At  least  for  the  idealized  conditions  considered  here,  the  New  Field 
Equations  describe  channel  tracking  about  as  well  as  the  complete  Frozen 
Field  Equations,  provided  the  axial  displacement  current  Is  retained  in  the 
New  Field  Equations.  Omitting  the  displacement  current  changes  tracking 
forces  qualitatively.  Physically,  oscillations  occurring  for  X  >  1  and 
early  time  magnetic  field  effects  are  lost.  Mathematically,  the  complex 
conjugate  pair  of  roots  are  replace  by  a  single  real  root  in  the  plasma 
dielectric  function.  (This  loss  of  a  mode  can  be  seen  more  simply  by 
deriving  the  dispersion  relation  for  normal  modes  of  the  New  Field  Equa¬ 
tions,  with  and  without  the  axial  displacement  current,  in  a  homogeneous 
resistive  medium.)  At  a  minimum,  the  axial  displacement  current  should  be 
kept  In  any  New  Field  Equation  calculation  of  tracking.  It  seems  prudent, 
more  generally,  to  retain  the  axial  displacement  current  whenever  the  con¬ 
dition  oa  >  2  is  not  well  satisfied. 
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